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ABSTRACT 

The observed slope at the high-mass end of the initial mass function (IMF) displays a remarkable 
universality in a wide variety of physical environments. We predict that competitive accretion, the 
ongoing accretion of gas from a common reservoir by a collection of protostellar cores, can provide a 
natural explanation for such a universal slope in star forming regions with metallicities Z > 10~^ Zq. 
In our discussion, we point out that competitive accretion will occur whenever a gaseous region 
has multiple Jeans masses of material and contains large-scale motions that are controlled by the 
gravitational potential well. We describe how and when these conditions can be reached during the 
chemical enrichment of the Universe, showing that they can occur for a wide range of metallicities 
and environmental conditions. We also discuss the ability of other physical processes to limit the 
effects of further accretion onto protostellar cores. Current theoretical and numerical studies show 
that competitive accretion is robust against disrupting effects - such as feedback from young stars, 
supersonic turbulence and magnetic fields - in all but the most extreme cases. 



Subject headings: stars: formation 
of state 



stars: mass function - early universe ~ hydrodynamics - equation 



1. INTRODUCTION 

The origin of the stellar initial mass function (IMF) 
remains one of the central topics of research in star for- 
mation theory and observation. The IMF describes the 
distribution of masses with which stars are born, and typ- 
ically takes a power-law form at the high mass end (> 0.3 
Mq) o i dN (X m~°'dm , with a taking a value of around 
-2.35 (|Salpeteij|1955h . The distribution of masses then 
turns over at lower masses, and is typica lly f itted with 
eithe r a series of broken power-laws ( e.g. iKro upa 2001,, 
|2002[) . or with a log- normal function (|Chabrierll2003l) . 

While there may be some observational debate re- 
garding details around the peak and turnover in the 
IMF (e.g. lElmegreen et allboOSD . the high-mass end 
is typically well described by the Salpeter slope, for 
a wide variety of environ mental co nditions. For ex- 
ampl e. Orion (|Hillenb ran d fc Carpen ter 200^ and p 
Oph (iLuhman fc RiekdiT998[ ) both exhibit Salpeter- like 
IMFs above ~ 0.5 M©, despite the fact that the Orion 
Nebula cluster contains a quartet of ionizing sources 
at its center, while p Oph is a relatively quiescent re- 
gion of star formation. Even in more extreme evi- 
ronments, such as R138 in Doradus (jMassev fc Hunted 
Il998f) - which also lies outside our own Galaxy in the 
somewhat lower metallicity Large Magellanic Cloud - 
the high-mass end of the IMF also appears to follow 
a Salpeter slope. Only for the globular clusters, the 
relics of extremely low-metallicity star formation, does 
the description of the high-mass end of the IMF start 
to become less clear, since the oldest clusters will have 
lost all stars > 0.8 Mq due to stellar evolution. In this 
case, we must re ly on more indirect measurements of the 
IMF d Parescc fc de Marchi"2000'; 'de Ma rchi et all 120001 : 
IChabrierii2003 : ide Grijs fc Parmentier.2007f ). with some 



studies suggesting a s omewhat flatter than Salpeter slope 
for masses > 0.5 M^dde Marchi et al.lll999t iKoch et all 
[200llde Marchi fc Pulonell2007D . 

The Salpeter slope is significant for two reasons. First, 
the fact that the power is negative shows that low-mass 
stars are significantly more abundant than their high 
mass counterparts. Second, the fact that the slope is 
steeper than -2 shows that most of the mass from the 
star formation process goes into the low-mass stars. Our 
current theoretical understanding of the first stars in the 
Universe - the so-called Population HI (or Pop. Ill) stars 
- suggest that this was not always the case, as these first 
stars are e xpected to have had very large masses (see e.g. 
iBromm fc Larson 2004; Glover 20Q5.). However, solar- 
mass stars with metallicities as low as [F e/H] = —5.9 
have been observed in the Galactic halo (|Frebel et al.l 
l2005l dOOl), and although they cannot by themselves 
constrain the full IMF at these very low metallicities, 
their existence is suggestive of an early transition from 
an IMF dominated by massive stars to one dominated by 
low-mass stars. 

For present-day star formation, there are two contrast- 
ing theories for the origin of the IMF. The first sug- 
gests that the vast majority of the final stellar mass 
comes from a local reservoir which envelopes the proto- 
star, an d is typically referred to as 'isolated' star for- 
mation (IShu. Li. fc Allenl l2004f ) or monolithic collapse 
(|McKee fc Tanll2003l ^. The mass contained in this reser- 
voir is set by some process that breaks the cloud up 
into dense structures, s uch as gravitational fragmenta- 
tion (e.g. iLarso nl '19731), turbulent fragmentation (e.g . 
iPadoan fc NOTc lhmd 2002: iHennebelle fc Chabriei<2008h 
or a comb ination of both (Klessen fc Burker tl 120001 : 
iKlessen et al. 2000; Klessen fc Burkcrt 20QI|). These iso- 
lated reservoirs then feed their protostar (or small N 
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group) without further influence from the surrounding 
cloud. 

Com petitive accretion (jZinneckeil Il982t iBonnell et al.l 
l2001aira ) presents an alternative picture, in which the 
high-mass slope of the IMF is the result of continued ac- 
cretion over and above the initial local collapse of the pre- 
stellar core. This model makes us e of the fact that th e 
majority of stars form in clusters (iLada fc Ladal [20031 ) . 
with the protocluster 'clump' providing the reservoir for 
on-going accretion. The dynamics of the gas and proto- 
stars, that naturally arise due to the gravitational poten- 
tial of the protocluster core, provide a distribution in the 
accretion rates that results in a distribution of protostel- 
lar masses. Stars which achieve masses larger than the 
typical Jeans mass in the cloud, which are those following 
the Salpeter-like section of the IMF, do so by accreting 
gas fro m a region which is larger tha n the average Jeans 
length (|Bonnell. Vine, fc Batefl2004[ ). 

In previous studies of star formation occurring in low 
metallicity gas at high redshift, attention has focussed 
primarily on the isolated star formation model. Typi- 
cally, these studies have sought to determine the charac- 
teristic mass scale for fragmentation of the gas through 
the consideration of its chemical and thermal proper- 
ties, either within the framework of a simplified dy- 
namical model (e.g. .Omu kai 2000; Omukai et al. 2005) 
or wi thin a full cosmological hydrodynamica l simulation 
(e.g. iSmith fc Sigurdssonl[2007l : iSmith et al.l[2009.1 . How 
this characteristic mass goes on to describe the type of 
IMF that we see today, one which covers a wide range in 
mass, has yet to be addressed. 

In this paper, we focus our attention on the competi- 
tive accretion model for the IMF, suggesting that it can 
provide a simple mechanism for a universal high-mass 
slope, over a wide range in metallicities. First, we give 
an overview of competitive accretion in Sj2l stressing that 
it is fairly insensitive to the conditions in the young clus- 
ter (such as initial density or velocity profiles). We also 
discuss the conditions required in order for competitive 
accretion to dominate over the isolated fragmentation 
models. In ^ we briefly discuss the observational ev- 
idence from local star-forming regions that favours the 
competitive accretion model. Based on the expected 
thermodynamic behavior of the gas, we show in 21 that 
the conditions for competitive accretion can be achieved 
readily in low-metallicity gas, suggesting that a Salpeter- 
like slope at the high-mass end of the IMF is a universal 
outcome of clustered star formation. In Sj5l we highlight 
a number of processes which could lessen the effects of 
competitive accretion, but argue that it is unlikely that 
these would be able to prevent the phenomenon entirely, 
provided star formation occurs in boimd clusters. 

2. AN OVERVIEW OF COMPETITIVE ACCRETION 

Competitive accretion is the process whereby individ- 
ual stars, contained in a bound group or cluster, compete 
to accrete from a communal reservoir of gas. This process 
requires that there exists a significant, distributed com- 
ponent of gas and that the gas and stars are mutually 
bound in the overall gravitational potential. Lastly, the 
gas needs to have some freedom to move in response to 
the local and global potential of the system. This last is 
unavoidable in a bound system where star formation oc- 
curs and any impediment such as turbulence or magnetic 



fields can only act to slow the accretion. (A discussion 
of such processes is given in 

In order for competitive accretion to generate an IMF, 
the accretion rates must necessarily vary considerably 
from star to star. This variation in the accretion rate 
occurs in the picture outlined above, since the accretion 
rate depends on a combination of the star's mass, the 
local gas density, the relative velocity of the gas and the 
star, and the position of the star in the gravitational po- 
tential of the full system. Bonnell et al. (2001a) pointed 
out that there are two regimes of accretion, depending 
on the relative gas-star velocity and the tidal field of the 
full cluster potential (i.e. that from the sum of the gas 
and stars). In the regime where the gas and stars are 
following similar trajectories, due to motion dictated by 
the cluster potential, then the cluster tidal field is the 
limiting effect and determines the accretion rate. This 
tends to be the case when protostars and the surround- 
ing gas are both infalling toward the cluster center, which 
typically occurs during the fragmen tation of th e gas, and 
thus the formation of new stars (e.g. IBonnell "et a l. 20081. 
If the stars are virialised or otherwise have velocities un- 
correlated to the gas, then the ability of a star to capture 
passing gas, known as Bondi-Hoyle accretion, dominates 
the accretion process. Finally, tidal accretion also de- 
scribes the mass evolution of the massive stars located at 
the centre of the cluster, as they are more or less static 
and feed off the gas that gets funneled to the center of 
the cluster, as it falls down the potential well. 

A forming protocluster is envisioned to start out in the 
tidal accretion regime, with newly formed fragments and 
gas having low relative velocities. The protostellar veloc- 
ities may also be low relative to the centre of mass of the 
system at this point. The accretion rate is a measure of 
a star's ability to attract gas to itself. The stars need not 
move through the system, sweeping up large volumes of 
gas, in order to have significant accretion rates. Instead, 
the gas responds to both the local and cluster potential 
and moves accordingly. Bonnell et al. (2001b) showed 
that for clusters with a p cx mass distribution, as is 
expected under gravitational collapse and as observed in 
stellar clusters, tidal accretion results in a fairly shallow 
IMF, with dN/dm oc m~^-^. 

Bondi-Hoyle accretion occurs after the stars have viri- 
alised, for example via interactions in the cluster center, 
generating large velocities relative to that of the infalling 
gas. A star's ability to accrete is then set by the star's 
mass and velocity in addition to the gas density 

m (X p . (1) 

Zinnecker (1982) first showed how Bondi-Hoyle accretion 
in a uniform density gas generates a Salpeter-like IMF 
dN/ dm (X m~^. For more realistic clusters with a density 
profile and an initial mass segregation, the resulting IMF 
is steeper with dN/dm cx m~^-^ (Bonnell et al. 2001b). 

Gas falling into a bound system's gravitational poten- 
tial has to follow the local potential gradient. This gas 
can only settle in local potential minima, which are cre- 
ated by the individual stars. Even in the cluster centre, 
the local gradient will be defined by the stars nearby and 
hence the gas will settle into an individual star's potential 
well. Gas cannot gather in regions where the gradient is 
not zero and hence cannot form new cores from unbound 
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gas. Ongoing star formation is limited to infalling cores 
that are already bound as they enter the cluster poten- 
tial, as otherwise t he tidal fields of the system will shred 
any unbound gas (jBonnell et al.|[2008[ ). 

There are a number of aspects of competitive accretion 
which often cause confusion. First, the stars do not need 
to move in order to participate in competitive accretion, 
as in this case it is the tidal radii of the stars inside 
the cluster which determine the accretion rates. Second, 
most stars in a cluster actually lose the competition and 
remain relatively low-mass stars (Bonnell & Bate 2006). 
Finally, the most massive star, forming in the centre of 
the cluster, accretes more through tidal accretion than 
Bondi-Hoyle accretion. Instead of moving at high rela- 
tive velocity, the most massive star remains in the centre 
of the potential - which it helps to define - and the gas 
flows down to it. 

The resulting IMFs are robust in terms of the cluster 
properties. A collapsing cloud which is not completely 
cold always develops a p oc density profile which 
generates a shallow IMF for the lower mass stars. The 
high-mass Salpeter-like IMF depends primarily on the 
nature of the accretion rate and, as discussed above, 
generates an IMF which varies from a = —2 to a = —2.5, 
depending on the density profile and any initial mass 
segregation. 

2.1. Fragmentation and timescales 

For competitive accretion to occu r requires tha t the 
cloud is able to promptly fragment (|Pringlel[l989D . and 
that the gas can accrete on timescales shorter than the 
gas removal timescale due to feedback, 

'n^orm — ^acc — '^emerge- (2) 

The first condition requires that the timescale of col- 
lapse of the individual fragment s be shorter than th e 
collapse timescale of the system (jBonnell et al.|[2b01al ). 
This will occur as long as the system contains many 
Jeans masses of gas and there is significant non-linear 
structure to act as seeds for the fragmentation. The 
important timescale here is the free-fall time for the 
system, tg, defined by tg = (37r/32G(p)voi)^^^, where 
(p)voi is defined to be the volume-av eraged density in 
the protocluster ([Bonnell et al.l[2b01a| ). This is also the 
timescale for stars to attain virialised velocities as they 
react to the cluster potential, and thus enter the Bondi- 
Hoyle accretion regime. Even if the cloud contains tur- 
bulent motions, the turbulent crossing time (tci) must 
be comparable to te, provided that the cloud is gravi- 
tational bound, since tg/tcr = (7ri?k)/(-\/T0 |-Ep|), where 
ii'k and \Ep\ are (respectively) the kinetic and gravita- 
tional eneries, for a uniform density spherical cloud. As 
such, turbulent motions typ ically do not alter the accre- 
tion within the cluster (e.g. lKlessenll2001l : iBonnell et al.l 
[200llSchmeia fc Klessenl l200l" note also that a more de- 
tailed discussion of the role played by turbulence is given 

m- 

The second condition is simply that the cluster of pro- 
tostars must be able to accrete some useful fraction of 
the background reservoir before t he reservoir is r emoved 
via feedback processes. Again, B onneU et al.l (|2001bD 
have shown that this is comparable to the mean free-fall 
timescale for the cluster as a whole. Naturally, additional 
processes may 'heat' the gas in the cluster's potential 



well, to the point where tg is no longer an accurate de- 
scription of how the ambient gas is accreted. However, 
current numerical simulations seem to suggest that this 
is not the case in young clusters, and a fuller discussion 
of these processes is given in The fact that tg is the 
governing timescale for both these conditions is part of 
the reason why the competitive accretion process is so 
robust once the general conditions are met, since it is 
not a strong function of the density. 

The fragmentation depends crucially on the thermody- 
namics of the gas , which determines the number of Jeans 
masses available ( LarsonllToSSL 120051 : iJappsen et al.ll2005l : 
IClark et al.l 12005 ). If the gas remains isothermal such 
that the cooling balances any heating, then even though 



the local Jeans mass, Mj cx 



1/2 



decreases during 



collapse, the number of Jeans masses available for frag- 
mentation remains constant. This is a consequence of the 
non-homologous nature of the collapse, which typically 
generates a density profile p cx r^^, leaving a central core 
of mass Mcorc oc p~^^^ where fragmentation can occur. 
In the p (X envelope, tidal forces will stop any frag- 
mentation. The major result is that the number of frag- 
ments produced is given by the number of Jeans masses 
available in the initial condit ions, with small adjustments 
due to the local geometry ( Larson 1985; Bastio n et al.l 
19911: iBur kert fc Bodenheimer 1996: ,Klessen fc BurkertI 
2000, 2001). 

If the initial conditions contain one Jeans mass, then 
the cloud needs to cool during collapse in order t o 
fragment (jRees fc Ostrikeil [1971 ^ lArcoragi et al.lll991h . 
whereas even if the system contains many Jeans masses, 
fragmentation will be increasingly suppressed if the cloud 
heats up with increasing density. Thus, fragmentation 
will predominantly occur in or follow from evolutionary 
regions where the temperature decreases with increasing 
density. Once fragmentation does occur, if it produces 
bound systems which can accrete, then competitive ac- 
cretion will necessarily follow and produce a ubiquitous 
shape to the IMF. 

It has also been demonstrated that competitive accre- 
tion is insensitive to the initial mass distribution of the 
protostars (or protostellar cores) that come out of the 
fragmentation process (Bonnell et al. 1997). Motions 
within the cluster, along with further accretion, tend to 
re-organize the protostars such that the resulting mass 
function produced in the cluster is the same. As such, 
competitive accretion will produce a Salpeter slope at 
the high-mass end that is robust against the details of 
how the gas is fragmenting. 

In summary, we suggest that any region of collaps- 
ing, fragmenting, gas, will undergo a competitive, Bondi- 
Holye-style accretion process. The small degree of varia- 
tion in the high-mass IMFs that this process produces is 
what we believe to be responsible for the near universal 
IMFs observed in all resolved stellar populations. In the 
following section we review the conditions in local star- 
forming regions, showing how they exhibit the properties 
that have been discussed here. In 21 we go on to discuss 
how these competitive-accretion-friendly conditions can 
be achieved for a wide range of metallicities. 



3. OBSERVATIONAL PICTURE FROM STAR-FORMING 
REGIONS 
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Given the picture outlined in ^ it is clear that 
for competitive accretion to be important in real star- 
forming regions, they must contain multiple Jeans masses 
and contain gas whose properties - such as the pri- 
mary fragmentation mode and motions - are described 
by timescales that are shorter than the time for the clus- 
ter to emerge from the cloud. Indeed, proto-clusters that 
are globally collapsing would provide a particularly fer- 
tile environment for competitive accretion. 

The results from recent observational surveys 
of ne arby star-formin g regi on s favor th is pic- 
ture (IWilking fc Ladal 119831: iMotte et all 119981: 




Testi fc SM-gentI 119981: iJohnstone et al.| 
iNutter fc Ward-ThompsonI l2007t 



2mi i2oniL 

lAndre et al.l 

Typically, the clouds contain upwards of 500 
Mq and the temperatures and densities are in the range 
of 10-20 K and IC^-IO^ cm~^, respectively. Such low 
temperatures are to be expected since the background 
interstellar UV is typically highly attenuated in regions 
of s tar formation due to their hig h column density 
Ce.g. IStamatellos fc Whitwor th 2003; lYoung et al1l2004l : 
iShirlev et al.ll2005D . The observed Jeans mass in the 
clouds is therefore < 2M0, and the corresponding Jeans 
length is < 10,000 AU. Not only do these clouds there- 
fore contain multiple Jeans masses, but the distance 
between pre-stellar cores is also observed to be similar to 
the Jeans length in the regions of active star formation. 
This suggests that the clouds are free to undergo an 
initial Jeans-like fragmentation. The free-fall time of 
the gas in these star forming regions is also quite short, 
of the order 10^ yr or less, much shorter than the 4 - 
5 million years it takes for the proto-stellar clusters to 
throw off their natal envelopes (jLeisawitz et al.l[l989t ). 

There is also evidence that these very young regions are 
not only multi ply Jeans unstable, but are also kinetically 
very cold. An dre et al.l ()2007| ) show that the virial mass 
derived from the N2H+ line emission in the L1688 cloud 
in p Oph is an order of magnitude small er than the value 
deriv ed from the C^^O observations (jWilking fc Ladal 
119831 ). Similar results have been found for other young 
regions such as NGC 2264 (Pcrctto ct al. 2006, 2007) and 
NGC 1333 (| Walsh et all 120061 120071 ). Such kinetically 
cold initial conditions are encouraging for the competi- 
tive accretion scenario, since they suggest that the ob- 
served pre-stellar cores will be brought together by the 
large-scale collapse of the region. 

4. THE DEPENDENCY ON METALLICITY 

Although the observations described in ^suggest that 
the gas in local star-forming regions is multiply Jeans 
unstable and hence in an ideal state for competitive ac- 
cretion to occur, we have no such direct probe of the 
state of the star-forming gas in lower metallicity sys- 
tems. However, from a theoretical viewpoint, it does not 
seem difficult to arrange for the gas to achieve the condi- 
tions required for competitive accretion. For competitive 
accretion to occur, we require multiple Jeans masses of 
gas, plus enough density substructure to act as seeds for 
fragmentation. If we start with ~ 1 Jeans mass of gas, 
then, as explained in ^2.11 in order to produce multi- 
ple Jeans mass, the gas temperature must decrease dur- 
ing the collapse, i.e. the gas must cool. Moreover, if we 
are to preserve our density substructure for long enough 
for it to cause fragmentation, then we must be able to 



cool quickly, within less than a dynamical time, since we 
do not expect gravitationally unbound substructure to 
persist for much longer than a dynamical time. For a 
gravitationally bound protocluster, the dynamical time 
tdyn is of the same order as the free-fall collapse time 
of the protocluster, tg, and so our requirement for com- 
petitive accretion is the same as the classical require- 
ment for efficient gravi tational fragmentation, tcooi < ts 
(iRees fc OstrikeiillQTTh . 

In local star-forming regions, we know that the cooling 
time is very short, and so this condition is easily sat- 
isfied. Can the same thing be said for low-metallicity 
gas? In the limiting case of zero metallicity, we know 
from the many detailed simulations that have been per- 
formed that this is generally not the case. Owing to the 
relative ineffectiveness of H2 as a coolant, the gas cools 
slowly, and achieves at best a cooling time that is compa- 
rable to its dynamical time, b ut that is not signifi cantly 
shorter (see e.g. Figure 6 in O'Sh ea fc Normaiil [2007). 
We therefore would not expect competitive accretion to 
occur in this case, and indeed numerical simulations 
that start with proper cosmological initial conditions 
find little evidence for widespread fr agmentation (e.g. 
lAbel et al.l [20021 : lYoshida et al.ll2006f ). although binary 
formation may occur in some cases (Turk et al., 2009, in 
prep.). However, the addition of even a small quantity 
of metals to the gas can dramatically alter this picture. 
The detailed chemical and thermodynamic evolution of 
low -metallicity gas has been the focus of several stud- 
ies fSchne ider et ani2002l[2006l : IOmukai et al.ll2005f ). and 
the general features of th eir results are broad ly similar 
(see, in pa rticular, Fig. 1 of lOmukai et al.ir2005l and Fig. 3 
of Schneider et al.ll2006l ). These studies have shown that 
there are two main phases during which metal-enriched 
gas can cool significantly, and on a timescale faster that 
the free-fall time. One occurs at low densities, and is 
caused by a combination of rotational and vibrational 
line cooling from H2 and HD and atomic fine structure 
cooling from carbon, oxygen, silicon and iron. The other 
occurs at higher densities, and is caused by dust cooling. 
These two regimes are illustrated schematically in Fig. 
1, and arc discussed in turn in the following sections. 

4.1. Low density cooling from line emission 

A number of stu dies, beginning with the work of 
iBromm et al.l (|2001|) . have examined the cooling of low 
metallicity gas caused by atomic a nd molecular line emis- 
sion at densities n < 1 0^ cm~'^ (iBromin fc Loebl[2?M 
Santoro fc Sluill [20061: iJappsen et al.l I2007L l2009allb: 
Smith fc SigurdssonI 120071: ISmith eFal] I2009D . These 



studies have shown that fine structure line emission from 
various metals, notably carbon and oxygen, is able to 
efficiently cool low density, metal-enriched gas, provided 
that the degree of enrichment is large enough. A com- 
mon way to quantify the required degree of enrichment is 
in terms of a critical metallicity Zcrit- This is, loosely, the 
lowest metallicity for which metal line cooling becomes 
an effective coolant at some point during the gravita- 
tional collapse of the gas. For instance, if we ask what 
metallicity is required in order for metal line cooling to 
balance compressional heating, we obtain an answer than 
depends on both the temperature and density of the gas, 
since the line cooling rates and the compressional heat- 
ing rate are functions of both of these variables. How- 
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Fig. 1. — We illustrate the expected equilibrium temperatures of the gas as a function of density for several levels of metal enrichment. Note 
that these are purely representative of the results from detailed studies. The solar metallicity line, for n < lO^'^cm"'^, is an approximation of 
that proposed by Larson ( 20051) . combined with the expected shift to an effective polytropic gamma of 1.4 at around n fs; lO^^cm"'^ once the 
gas become optically thick (e.g. IMasunaga. Mivama. &: Inutsukalll998l ). and then the subsequent cooling provided wrhen the temperature is 
high enough to dissociate H2. The other three curves for the lower metallicity and primordial gas arc taken to follow (roughly) the results 
from the study by Omukai et al (2005), in particular their Fig. 1 (but see also .Schneider ct al.i.2002i . i2006i) . All the metal enriched curves 
show two clear regimes: one in which the cooling is dominated by the fine-structure line-emission (for n < lO^cm"'^, depending on the 
metallicity), and a second in which the cooling is dominated by the dust. For lower metallicity gas the dominant formation process for 
H2 is via three-body gas phase reactions, in contrast to the solar metallicity case in which the dust fraction is large enough for the grain 
surface reactions to dominate. The energy released during the three-body reactions offsets the cooling that the H2 provides, and as such 
the gas is unable to rid itself of the compressional heating caused by the gravitational collapse, resulting in the temperature rise between 
the two cooling regimes. Fragmentation has been shown to occur in the regimes where the temperature decreases with increasing density. 



ever, the gas does not evolve arbitrarily in temperature- 
density space, and along any realistic trajectory, Zcrit al- 
ways has a well-defined minimu m that falls in the range 
10-''Z(T, < Z rrH. < IQ-^-^ (iSantoro fc Shulll 120061 : 



iFrebel et al.l[2TOl [siuig[200l 

Numerical simulations have confirmed that gas 
with a metallicity Z > Zcrit can cool efficiently and 
quickly at low densities until the temperature hits 
the floor set by the CMB tem perature, Tqmb = 
2.726(f + z) dB romm et al.' '200f'; ^Smith fc Sigurdssonl 
l2007HSmith"et al. 2009; Jappscn ct al. 2009a, l^. The gas 
is therefore able to increase the number of Jeans masses 
available for gravitational fragmentation, while preserv- 
ing the substructure required to seed that fragmentation. 
We therefore potentially have ideal conditions for com- 
petitive accretion. However, in practice, there is a com- 
plication: if the initial gas mass in the cooling region was 
significantly less than a Jeans mass (as may be the case 
if the gas is initially confined by a gravitational potential 
dominated by dark matter, rather than by its own self- 
gravity), then the gas may not be able to cool to a low 
enough temperature to produce multiple Jeans masses, 
and hence may not fragment. In local star-forming re- 
gions, this is not generally a concern, as the gas tem- 
perature can reach very small values. However, at high 
redshift, the temperature floor set by the CMB is much 
higher, and hence this is far more of a concern. 

The various numerical studies that have been per- 
formed of metal-enriched gas evolving in early proto- 



galaxies tend to support this concern. Simulations 
that start with idealized in itial conditions (Bromm et al.l 
I2OOII : iJappsen et al.|[2009a[) produce large masses of cool- 
ing gas that can easily become multiply Jeans unstable 
and that readily fragment. On the other hand, simula- 
tions that follow the build-up of a metal-enriched pro- 
togalaxy from more realistic cosmological initial condi- 
tions tend_to_2rodjj£e_amt^^ smaller mass of cooling 
gas ()Smith fc Sigurdsso^ 120071 : [S^ith et all (2009). In 
these simulations, fragmentation occurs only for metal- 
licities in the range 10"^-^ Z© < Z < 10"^ '^ Zq. At lower 
metallicities, metal-line cooling never becomes impor- 
tant, and hence icooi remains comparable to or larger 
than tff throughout. On the other hand, at metallicities 
Z > 10"^-^ Zq the gas cools too quickly: it reaches the 
CMB temperature floor at an early stage in the collapse, 
while its Jeans mass is still comparable to the amount of 
mass contained within the cooling region. Thereafter, it 
evolves isothermally at Tcmb, and hence cannot further 
increase the number of Jeans masses ava ilable for frag- 
mentat ion. Similar results were found bv IJappsen et al.] 
i2009b) in a study of Z = 10 Z© gas cooling from an 
initially ionized state: the enhanced ionization promoted 
rapid cooling, and allowed the gas to reach the CMB 
temperature floor while still at low density, the number 
of Jeans masses available for fragmentation at the tran- 
sition to isothermal evolution was small, and so the gas 
did not fragment. 
The important role played by the CMB in regulating 
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low density fragmentation was furt her highlight e d by a 
numerical experiment performed by lSmith et alJ (|2009r ). 
They ran a simulation with Z = 10~^ Zq in which the 
thermodynamical effects of the CMB were not included 
(i.e. the CMB temperature was artificially set to zero). In 
this case, the gas was able to cool down to around 20 K 
(compared to 70 K in the case with the CMB, as the 
simulated collapse occurred at a redshift z = 25), allow- 
ing the number of Jeans masses of gas contained within 
the cooling region to be increased by a further factor of 
ten. This resulted in fragmentation, as one would expect, 
with a larger number of fragments forming in this case 
than in any of their other simulations. 

The low density, line-dominated cooling regime does 
not extend far beyond n — 10® cm~^. Above this den- 
sity, all of the major fine structure coolants have reached 
local thermodynamic equilibrium (LTE), and hence the 
cooling time due to these processes tends to a constant 
value, while the free-fall time continues to decrease with 
increasing density. We therefore soon reach a regime in 
which icooi > iff, at which point compressional heating 
starts to warm up the gas. In addition, H2 formation 
heating also begin s to make a significa nt contribution at 
around this point (jOmukai et al.ll2065[ ). The cloud there- 
fore enters a 'loitering' phase, allowing material to collect 
at these densities. 

4.2. High density cooling from dust continuum emission 

The loitering phase that begins once metal-line cooling 
can no longer cool the gas effectively continues until the 
onset of the second phase of efficient cooling. This occurs 
once the timescale for thermal energy transfer from the 
gas to the dust grains becomes small. To determine when 
this occurs, we need to know how the properties of the 
dust - specifically, the grain size distribution and the 
ratio of dust mass to gas-phase metal mass - compare 
to those of the dust present in local star- forming regions. 
The rate at which energy is transfered from gas to dust is 
proportional to the total available surface area of grains 
per unit volume of gas, Atot- We can write the ratio of 
Atot in gas with metallicity Z to Atot in solar metallicity 
gas as 

Aot(Z) _ „ ^ 

A (rj \ ^ /dust „ ■ \'J) 

If we assume that the dust in low-metallicity gas is sim- 
ply a l ess abundant version of local dust ( i.e. that it re- 
tains a lMathis. Rumpl. fc Nordsieckl (|1977l ) size distribu- 
tion and that the ratio of dust mass to gas-phase metal 
mass is the same as in local gas), then we would simply 
have /dust = 1, and in this case the onset of the second 
phase of efficient cooling begins at a density of ro ughly 
ricr 109(Z/I0-3 Zo)-^ cm-3 (lOmukai et aLlIM^ . In 
the more general case, where the dust need not have 
these properties, we can instead write the required den- 
sity as ricr ~ 10Vdust(Z/10"^ Zq)-^ cm'^. The actual 
value of /dust in high redshift, low metallicity systems is 
not known, and is an active field of research. The large 
quantities of dust ob served in the host galaxies of z ~ 6 
quasar ho st galaxies (iBertoldi et aT]|2003t[Maiolino et al.l 
[2004; Wa ng et al.l[200l " suggest that some mechanism 
for efhciently producing large quantities of dust must be 
active in the early universe. On the other hand, the 
latest theoretical models of dust production in super- 



novae - formerly the leading candidates for producing 
the observed high-redshift dust - suggest that much of 
the dust that is initially formed in the supernova rem- 
nant is destroyed by the passage of the reverse shock, and 
that little survives to escape into the in terstellar medium 
(|Nozawa et all 120071: iBianchi fc Schneider .2007.) . It is 
therefore quite plausible that /dust "C 1 in low metallic- 
ity gas, but equally it is also plausible that /dust ~ 1- 

Once the transfer of energy from the gas to the grains 
becomes rapid, then the gas quickly cools down to 
the dust temperature T^. This sudden cooling occurs 
throughout the large quantity of gas accumulated during 
the loitering phase, making it particularly easy in this 
case to produce a highly Jeans unstable gas distribution. 
An example of t his mechanism in action can be found in 
iClark. Glover, fc Kless e n (i,2008i). whe r e the temperature- 
density relations from FOmukai eTaH (I2OOI are approx- 
imated by a piecewise-barotropic equation of state. It 
was found in this study that vigorous fragmentation oc- 
curs during the dip in the equation of state caused by the 
sudden onset of efficient dust cooling, resulting in a clus- 
ter in which competitive accretion processes dominate. 
Indeed, the result ing mass function of the 'sink particles' 
()Bate et al.lll995[ ) was consistent with a Salpeter slope 
for metallicities as low as Z = 10^^ 7jq. 

4.3. Implications for the IMF 

As the preceding sections show, there are two regimes 
in low-metallicity gas in which competitive accretion 
could possibly operate: a low-density regime in which 
the cooling is dominated by atomic and molecular line 
emission, and a higher density regime which is dominated 
by dust cooling. In local star forming regions, these two 
regimes almost overlap: molecular line cooling begins to 
become ineffective at densities n ~ lO'' cm~'^ and tem- 
peratures T ~ lOK, owing to the large line opacities that 
have developed by this point, while dust cooling starts to 
becomes highl y effective at only a slightly higher density 
n - 10^ cm-3 (|Goldsmithll200ll) . 

In lower metallicity gas, these two regimes are more 
distinct. Their relative importance depends on two 
main factors: the metallicity of the gas, and the de- 
gree of heating provided by the cosmic microwave back- 
ground or some other background heating source. At 
the lowest metallicities, neither metal-line cooling nor 
dust cooling is able to cool the gas rapidly enough to al- 
low competitive accretion to occur, and the evolution 
of the gas takes place just as in the metal-free case 
(|Clark. Glover, fc Klessenl[200l . 

At higher metallicities, dust and/or metal-line cooling 
become effective. The required metallicity depends on 
the properties of the dust. Dust cooling becomes effec- 
tive at a metallicity Z ~ 10"^/^^^^ Z©, where /dust quan- 
tifies the total available surface area of dust grains, as 
described in ^ 34.21 above, whereas metal-line cooling be- 
comes effective for a critical gas-phase metallicity Zcrit = 
10~'^'^ Tjq. Therefore, if /dust ^ 0.03, metal-line cooling 
becomes effective at a lower metallicity than dust cool- 
ing. In this regime, the gas cools faster than it can dy- 
namically respond. However, it cannot cool indefinitely; 
eventually, it must reach a temperature at which further 
cooling is balanced by heating, either from compression, 
or from the effects of external sources of radiation or 
the CMB. The Jeans mass of the gas at the point where 
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this phase of efficient cooling finishes - due, for instance, 
to the gas temperature reaching the CMB floor - de- 
termines whether competitive accretion will be possible. 
It sets the characteristic mass scale for any protostellar 
cores that form, and since competitive accretion requires 
a collection of such cores, it requires that the mass of 
cooling gas available be much larger than the Jeans mass 
at this point. If this is the case, then the outcome will be 
an IMF with a power-law high-mass end, and a turnover 
corresponding to this characteristic mass scale. 

For example, if metal-line cooling is halted by the gas 
temperature reaching the CMB temperature, and this 
occurs at a gas density noM R, then this character istic 
mass scale can be written as (jElmegreen et al.]|2008[) 

'rj^ X 3/2 ^ 

J CMB \ / "•CMB 



50 K 



/ "■CMB \ 

V 104 cm-3 / 



-1/2 



Mr. 



(4) 



Typically, ncMB ~ 10** cm"'^ for metallicities close to 
Zcrit, and so at high redshift, we expect Mchar in these 
systems to be much larger th an a solar mass (c.f. the 
fragment mass di stributions in lJappsen et al.ir2009al and 
ISmith et al.ll2009l ). If the metallicity is much larger than 
Zcrit, however, so that cooling is very efficient and "cmb 
is small, then Mchar may be comparable to the total mass 
of cooling gas, in which case gravitational fragmenta- 
tion will be strongly suppressed and competitive accre- 
tion will not occur. On the other hand, at low redshifts, 
Mchar ~ 1 Mq or less, and so in this case, we would 
expect a standard Salpeter IMF to be produced at all 
metallicities greater than the critical value. 

If /dust^O.03, then dust cooling will become effective at 
a lower metallicity than metal-line cooling. In this case, 
rapid cooling, fragmentation and competitive accretion 
will take place at high densities, where the characteristic 
Jeans mass is low even if heating from the CMB keeps the 
dust temperature at a much higher level than is typical of 
local star-forming regions. Therefore, we would expect 
a standard Salpeter IMF extending down to relatively 
small masses to be produced in this scenario. 

Once both dust and metal-line cooling have become 
important, the outcome depends on just how effec- 
tive the metal-line cooling actually is. For metallicities 
Z < 0.01 Z0, the metal-line cooling phase is followed 
by a he ating phase prior to the onset of efficient dust 
cooling ( Omukai et al.l[2005[ ). allowing the conditions for 
competitive accretion to be achieved for a second time 
during the same collapse, only now occurring in gas with 
a higher density, and hence lower characteristic mass. In 
this case, we would again expect a standard IMF. On the 
other hand, if Z > 0.01 Zq, then the metal-line cooling 
will keep the gas close to the CMB temperature until it 
becomes optically thick. In this case, even if dust cool- 
ing becomes effective at high densities, it cannot cause 
a significant drop in temperature; the gas continues to 
evolve isothermally, and competitive accretion does not 
occur. We would therefore expect to obtain a power- 
law IMF with Mchar ^ IM0 at high redshifts, smoothly 
transitioning to a standard IMF at low redshifts as Tqmb 
decreases. 

Finally, it should be noted that we are not the first 
to suggest that there may be a lin k between the CMB 
temperature and Mch a r (see e.g. iLarsonl [19981 [2OO5I : 



IClarke fc Brommll20Q3l : lTumhnsonll2007h . However, pre- 
vious work on this issue has generally focussed only on 



fragmentation during the metal-line cooling regime, and 
has concentrated on determining Mchar, rather than on 
determining the shape of the IMF at M > Mchar- 

5. PREVENTING COMPETITIVE ACCRETION 

Given any cloud with multiple Jeans masses, in which 
the gas and the newly formed fragments are free to move 
in the potential, competitive accretion will occur. Essen- 
tially, the gas will always want to reduce its energy state 
by moving down the potential well. Preventing compet- 
itive accretion requires that the gas reservoir be held up 
against the gravitational force of the cluster, but, at the 
same time, some of the gas must be permitted to form 
cores, or else no star formation would occur. In other 
words the ambient gas would need to behave as if it is 
'hot', while the cores are 'cold'. In this section we dis- 
cuss some of the processes that may lessen the effects of 
competitive accretion in real clusters. 

5.1. Supersonic turbulence 

Although most of the studies to have looked at com- 
petitive accretion have used calculations in which the 
supporting supersonic turbulent motions have been free 
to decay by dissipating energy in shocks, several comple- 
mentary studies have looked at t he cas e of driven tur- 
bulence. In particular, Klcssen' ()2001h examined how 
the turbulent driving scale affects the fragmentation of 
self-gravitating isothermal gas in multiply Jeans unsta- 
ble clouds. It was found that when the turbulence was 
driven on scales larger than the Jeans length, the star for- 
mation occurred in clusters. In contrast, when the tur- 
bulence was driven on progressively smaller scales, star 
formation became a progressively more isolated event, 
with competitive accretion playing almost no role. In 
the clustered mode, the mass function was set by the 
competitive accretion which occurred within the clus- 
ter's potential, and was consistent with the field star 
IMF (see also iJappsen et ahl I2005D . In the later case 
of isolated star formation, competitive accretion played 
no role and the resulting mass function was considerably 
flatter than the observed distribution. A more general 
discussion of the role of s upersonic turbulence in sta r 
formation can be found i n iMac Low fc KlessenI (|2004l ): 
iBallesteros-Paredes et~all ()2007( ) 

That competitive accretion can work under the condi- 
tions of large-scale driving in isothermal gas is not sur- 
prising: the ram pressure is not isotropic and entrains 
mass, such that the flows gather Jeans unstable material 
together to form large potential wells. Also, the turbu- 
lence is free to decay inside these over-densities, as it is 
typically shielded from the outside driving by a bound- 
ary shock. Once inside the potential, the gas is free to 
move under its influence when it has lost enough of its 
kinetic energy in shocks. This loss of kinetic energ y oc- 
curs on the crossing time ()Mac Low et al.lll998al lbf) and, 
as discussed above, can help to promote fragmentation. 
The fact that turbulence is observed to be dominated by 
large -scale modes suppo rts this picture of cluster forma- 
tion (IMac Low fc Ossenk opi 2000; Osscnk opf et al.ll2001l : 
lOssenkopf fc Mac Lowll2002c ,Brunt„200a) . 

5.2. Small-scale sources of turbulence 

In real clouds it is still not clear whether or not su- 
personic turbulence is also constantly driven on smaller 
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scales within the cluster, in a way which could dilute 
the effects of competitive accretion - putting aside for 
the moment that it appears to result in the wrong IMF 
(Klessen 2001). Jets and outflows from young proto- 
stars perhaps provide the best candidate for halting 
subsequent accretion (e.g. iMatzneil l2007l ). but a full 
treatment of the prob lem has never been implem ented. 
iLi fc Nakan^ (|2006D and lNakamura fc Lil ()2007D have 
looked at the effect of outflows (winds and jets) on the 
star formation rate in small clusters, concluding that 
they help it to self-regulate at a fairly low value. However 
these calculations do not attempt to model any subse- 
quent accretion onto the 'star' particles after their birth. 
Their numerical set-up also employs periodic boundary 
conditions, such that the outflow energy c an never escape 
the cluster. In a different set-up, Banerjee et al.l ()2007f ) 
performed high-resolution simulations of single jets, con- 
cluding that they were unable to generate space-filling 
supersonic turbulence. Rather they found that the tur- 
bulence was only supersonically driven in the bow-shock 
region, while the rest of the cloud exhibited sub-sonic 
turbulence. 

Altogether, it is not clear from the current simulations 
that competitive accretion in such an environment would 
be halted, rather than just slowed down. The results 
of Dale & BonncD (2008), which focus on the impact of 
winds from high mass stars, suggest the latter: the accre- 
tion/fragmentation process is slowed, but the winds have 
little effect on the mass function produced by competitive 
accretion. They find that both filaments and disks will 
very effectively collimate the outflows, allowing accretion 
to continue along the filaments. 

Ideally, jets and outflows have a better chance of kinet- 
ically heating the cluster gas reservoir - and thus slowing 
down the competitive accretion process - when there are 
multiple sources present in the cloud, since they will de- 
posit the ir combined en e rgy ov er a wider volume. For 
example, iBaneriee et al.l (|2007[ ) point out that it would 
take ~ 100 jets to match the kinetic energy of a virialized 
cluster such as the Orion Nebula Cluster (ONC), assum- 
ing a perfect coupling between the jet and the cloud. If 
one assumes a constant rate at which stars are formed, 
tsF, and a constant lifetime of the outflow, tjot, then 
the number of outflow sources available at any one time 
is, iVjet ~ tsp/ijot- For a system like the ONC with 
~1000 stars, and taking the typical lifetime of a jet to 
be 10^ yr, this implies a tsF of less than 10® yr, faster 
than some time scales for the ONC in the literature (e.g. 
iTan et ahl I2006D . In practice, a better measure of tjet 
might be the time at which it leaves the main cluster- 
forming core, since it will have little influence after this 
point. In addition, it has been suggested th an the jet 
hfetinie in the OMCl/2 is only a few 10^ yr (|Aso et al.l 
120001 IWilhams et"alll2003f l. If jets are to be important, 
then these shorter outflow lifetimes would suggest an 
even shorter cluster formation time. Finally, early results 
from self-consistent simulations of jets/outflows - which 
include coupling between the measured accretion rate 
and the launch velocity - suggest that the cloud is not 
substan tially disrupted, and accretion is only marginally 
slowed (jBanerjee et al.l[2009t ). 

As one goes to lower metallicity, the picture becomes 
even more uncertain. On the one hand, the winds 
from lower metallicity stars are expected to be weaker 



than those we see from solar metallicity star formation 
(|Kudritzkill200l IVink fc de Koteill2005h . On the other 
hand, it is likely that the thermodynamic behavior at 
lower metallicities leads to a higher rate of star forma- 
tion, as defined by the rate at which new stars form 
(^sf), and so their combined effect may be greater due 
to a better volume coverage. This is especially true 
if the fragmentation is dominated by the dust-coupling 
regime, since it typically occurs at higher densities than 
the line cooling, and hence is described by a smaller free- 
fall time. However, for the winds to be effective, they 
still have to counteract the large-scale infall of the gas 
that is coming in from lower densities, and destroy the 
filamentary structure that it wil l posses from the large- 
scale gravitational instabilities ( Klessen fc Burkertl |2000l : 
IVdzq uez-Semadeni et al. 2007,; Bonnell et al. 20081). In 
contrast, if the fragmentation takes place primarily in 
the line-cooling regime, then we would have a similar 
situation as in solar-metallicity star formation. 

5.3. Radiative feedback from young stars 

To date, the majority of the simulations used to study 
competitive accretion have neglected the effects of the 
accretion luminosity that is released as the gas falls onto 
the protostellar surface. Given that this radiat ion can al- 
ter th e Jeans mass in the surrounding gas (e.g. lKrumhoI3 
|2006[ ). and thus change the fragmentation characteristics 
around the protostars, it is instructive to ask whether 
this process can prevent competitive accretion. 

In the picture described in ^ all stars start off with 
masses comparable to the mean Jeans mass in the cloud, 
and so competitive accretion begins a s low-mass st ar for- 
mation. The recent calculations by iBatd (|2009f ) show 
that the heating from this low-mass star formation phase 
is important on scales of 100-1000 AU, much shorter than 
either the mean Jeans length in the parent cloud or in 
nearby star-forming regions. While this cuts down the 
number of objec ts forming in disks (in comparison to 
iBate et aI]|2003D , the larger-scale fragmentation is unaf- 
fected, and so is the dynamical interaction of the stars 
in the potential. 

However as some of the stars start to accrete more 
mass, the accretion luminosity will start to rise, par- 
ticularly if the accretion rate is high. For example 
iKrumholz. Klein, fc McKed ()2007al ) show that for accre- 
tion rates between 10""* to 10"'^ Mq yr~^, the gas can be 
heated to as much as 50 K at distances of around 2000 
AU. By this point in the competitive accretion picture 
however, the majority of the new fragmentation is taking 
place in dense filaments which are squeezed as they fall 
into the potential (BonneU, Clark, & Bate 2008). Not 
only are these filaments well shielded by the radiation 
from the central accreting sources, but they become grav- 
itationally unstable at distan ces > 8000 AU, still well 
outside the 50 K gas. Further. iDale et al] ()2005l ) showed 
that these filaments are also shielded from any ionizing 
sources at the cluster center. 

At lower metallicities, and in the case where the only 
period of competitive accretion occurs in the regime dom- 
inated by atomic fine structure cooling, a similar argu- 
ment applies. Fragmentation in this case occurs at sim- 
ilar densities to those at which it occurs in present-day 
star- forming clouds, and since the characteristic mass is 
comparable to or larger than the present-day value, the 
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associated length scales will also be comparable to or 
larger than those in present-day clouds. 

If fragmentation occurs in the regime dominated by 
dust cooling, however, then the situation is somewhat 
more complicated. As one moves to lower metallicities, 
the coupling to the dust ~ and thus the fragmentation 
- occurs at progressively higher densities. The question 
is then whether the heating of the dust by the accre- 
tion luminosity is sufficient to remove the dip in the 
temperature-density curve (see Figure 1). As the dip 
is caused by the gas cooling to reach the dust tempera- 
ture, the dip can be removed if the heating due to the 
accretion luminosity raises Tdust above Tgas before the 
onset of efficient thermal coupling between gas and dust. 
How easy this is to bring about depends on the gas tem- 
perature before the onset of efficient thermal coupling. If 
the gas temperature is high, as will generally be the case 
in low metallicity systems, then the dust must also be 
heated to a high temperature. Since the rate at which 
the dust can re-radiate energy sc ales as a high power 
of the dust temperature (typically Adust c>c T^^^^ - r^ust 
for Tdust < 100 K), the energy input required in this 
case is very large, and is only likely to be achieved in 
gas very close to the accreting protostar. More plausi- 
bly, heating due to accretion may reduce the depth of the 
dip in the temperature-density curve slightly, thereby in- 
creasing the characteristic mass. However, if cooling and 
fragmentation are occuring at high densities, then the 
characteristic mass will remain relatively small. If the 
metallicity is higher, and the gas temperature is lower, 
then increasing Tdust to Tgas requires less energy input, 
but at the same time, the thermal coupling and frag- 
mentation will occur at lower densities, further from the 
protostar, and hence in a region where there is far less 
heating. 

5.4. Magnetic fields 

Magnetic fields are ubiquito usly observed i n in- 
terstellar gas on al l scales (jCrutcher et al.l 120031 : 
iHeiles fc Trolandll2005( ). However, their importance for 
star formation and for the morphology and evolution of 
molecular cloud cores remains controversial. A crucial 
parameter in this debate is the ratio between core mass 
and magnetic flirx. In supercritical cores, this ratio ex- 
ceeds a critical value and collapse can proceed. In sub- 
critical cores, magnetic fi elds provide stabihty (Spitze3 
Il978t iMouschoviasHlQQlal bl) at least if we assume ideal 
magnetohydrodynamics. Recent measurements of the 
Zeeman splitting of molecular lines in nearby cloud cores 
indicate mass-to-fiux ratios that lie above the critical 
value, in some cases only by a small margin but very 
often by factors of many if non-detection s are included 
(ICrutcheilll999HCrutcher fc Trolandll2000t iBourke et all 
120011: iCrutcher et al.ll2008l )! 

We now consider the conditions for competitive accre- 
tion, as discussed in Section 12. 1[ in the context of su- 
percritical and sub critical regimes. For subcritical cloud 
cores to form stars (jShu et al.lll987f l. ambipolar diffusion 
is required, i.e. the drift between neutral and charged 
particles. It causes a redistribution of the magnetic flux 
until the inner regions of the c ore become supercritical 
and go into dynamical c ollapse (Mouschovias 1976l ll979l : 
iMouschovias fc Paleologou 1981 ). This process was orig- 
inally thought to be slow, because in highly subcritical 



clouds the ambipolar diffusion timescale, tAD, is about 
10 times larger than the dynamical time, tg. Regardless 
of the exact value of the ambipolar diffusion timescale, 
in a magnetically dominated region the first two system 
timescales (rtorm, and race) are roughly equivalent to iAD- 
Their relative values remain unchanged and the first rela- 
tion in condition [5] still holds. However, the last relation 
may break down for tAD ^ iff , and feedback may remove 
the cluster gas before contracting cores have had time to 
interact or undergo further accretion from the common 
reservoir. 

However for cores close to the critical value, as is sug- 
gested by observations, both timescales are comparable. 
Numerical simulations furthermore indicate that the am- 
bipolar diffusion timescale becomes significantly shorter 
for turbulent velocities similar to the values observed 
in nearby star-forming region dFatuzzo fc Adams! [20021 : 
iHeitsch et all 120041: iLi fc Nakamural I2004D . In this case 
we fully recover both relations in Expression [21 and the 
star forming regions should be subject to competitive ac- 
cretion. If most cloud cores are magnetically supercriti- 
cal anyway, then the relevant timescale is the dynamical 
one, iff, and the discussion in Section [2. II is applicable. 

However, we note several points of caution. Although 
magnetic fields in general appear to be too weak to pre- 
vent gravitational collapse to occur, they still may influ- 
ence the star formation process in various ways. Recent 
numerical simulations have shown that even a weak mag- 
netic held can change the coupling bet ween stellar feed- 
back processes and their parent clouds (jNakamura fc Lil 
12007; Krumholzjjt a l. 20 0713 ), or slow down the over- 
all evolution ( Heitsch et al] |2"001[ ) by a factor of two 
or so. Although this increases Tacc, it will, under nor- 
mal conditions, remain smaller than Temerge- Com- 
petitive accretion still influences the evolution but at 
a lower level than in the purely hydrodynamic case. 
We also mention, that magnetic fields when considered 
in the ideal MHD approximation are able to reduce 
t he fragmentation of cloud cores on very small scales 
{Price & Bate"2007b', '2008'; 'Hcnncbclle_&iromang"200l; 
HcnncbcUc fc Tcvssicr 2008) and infl uence the properties 
of pro tostellar disks (P rice fc Bate! l2007al: [Mellon fc Q 
l2008al |bl). Due to magnetic breaking, disks tend to be 
smaller and less prone to gravitational instability. This 
efficiently suppresses binary formation in such models, 
and again points towards the necessity to consider non- 
ideal MHD (such as ambipolar diffusion or Ohmic dis- 
sipatio n) in order t o explain the observed binary distri- 
bution ()Ladall2006[ ). Our current analysis is only weakly 
affected by these uncertainties, because competitive ac- 
cretion acts on scales much larger than individual proto- 
stellar disks. 

6. SUMMARY 

Current studies of resolved stellar populations flnd the 
that high-mass end of the initial mass function (IMF) is 
well described by a Salpeter-like slope (Salpeter 1955). 
In this paper, we predict that whenever star formation 
occurs in bound stellar clusters, the slope of the mass 
function, above some characteristic mass, will be consis- 
tent with the Salpeter value. We suggest that this is due 
to the process of competitive accretion (e.g. Bonnell fc 
Bate 2006). 

In our discussion, we point out that competitive ac- 
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crction is robust against both the protocluster density 
profile, and the details of the fragmentation within the 
protocluster's gravitational potential. The only require- 
ment is that fragmentation can result in the formation 
of a bound group of protostars. 

We explore how fragmentation is likely to vary as a 
function of metallicity and show that for non-zero metal- 
licities, there exist the necessary regimes where gas cools 
while collapsing, such that a single Jeans mass cloud will 
evolve towards having numerous Jeans masses. Depend- 
ing on the metallicity of the gas, this cooling can be pro- 
vided by either line emission, occurring at 'low' densities 
(n < lO^cm^'^) or dust continuum emission, occurring at 
higher densities. At solar metallicities, the two regimes 
coincide. The result of this cooling is fragmentation of 
the collapsing gas cloud, and the formation of a bound 
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